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Abstract. We recover the Navier-Stokes equation as the incompressible limit 
of a stochastic lattice gas in which particles are allowed to jump over a meso- 
scopic scale. The result holds in any dimension assuming the existence of a 
smooth solution of the Navier-Stokes equation in a fixed time interval. The 
proof docs not use non-gradient methods or the multi-scale analysis due to the 
long range jumps. 



1. INTRODUCTION 

A major open problem in non-equilibrium statistical mechanics is the derivation 
of the hydrodynamical equations from microscopic Hamiltonian dynamics. The 
main difficulty in this project lies in the poor knowledge of the ergodic properties of 
such systems. To overcome this obstacle, deterministic Hamiltonian dynamics have 
been successfully replaced by interacting particle systems (cf. fJ* and references 
therein) . 

Following this approach, in the sequel of the development of the non-gradient 
method by Quastel and Varadhan [6], Esposito, Marra and Yau [H [2] derived 
the incompressible Navier-Stokes equation for stochastic lattice gases in dimension 
d>3. 

The main step of their proof relies on a sharp estimate of the spectral gap of 
the jump part of the generator of the process and on the characterization of the 
germs of the exact and closed forms in a Hilbert space of local functions. The 
characterization of the closed forms as the sum of exact forms and currents allows, 
through a multi-scale analysis, the decomposition of the current as a sum of a 
gradient part and a local function in the range of the generator. 

In this article we consider a stochastic lattice gas with long range jumps. The dy- 
namics is build in a way that the density and the momentum are the only conserved 
quantities. Choosing appropriately the size and the rates of the jumps, we are able 
to show that a small perturbation of a constant density and momentum profile 
evolves in a diffusive time scale as the solution of the incompressible Navier-Stokes 
equation. 

In contrast with [HE], the mesoscopic range of the jumps permits to consider 
perturbations around the constant profile of order N^'', for b small, where N is 
a scaling parameter proportional to the inverse of the distance between particles. 
This choice has two important consequences. On the one hand, in order to close the 
equation, one does not need to replace currents by averages of conserved quantities 
over macroscopic boxes, but only over mesoscopic cubes, whose size depend on the 
parameter b. In particular, there is no need to recur to the multi-scale analysis or 
to the closed and exact forms, simplifying considerably the proof. On the other 
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hand, choosing b small enough (6 < 1/2), one can avoid in dimension 1 and 2 
the Gaussian fluctuations around the hydrodynamic limit and prove a law of large 
numbers for the conserved quantities in this regime. We are thus able to derive 
the incompressible Navier-Stokes equation even in low dimension, where the usual 
approach is intrinsically impossible since it involves scales in which fluctuations 
appear. 

The main drawback of the approach presented is that it requires a bound on 
the spectral gap of the full dynamics restricted to finite cubes. The bound needs 
only to be polynomial in the volume of the cube, but the generator includes the 
collision part. This problem, already mentioned in [2J, is rather difficult in general. 
We prove such a bound in Section [6] for a specific choice of velocities. 

The model can be informally described as follows. Let V be a finite set of 
velocities in R'', invariant under reflections and exchange of coordinates. For each v 
in V, consider a long-range asymmetric exclusion process on Z'' whose mean drift is 
vN~^'^^''\ Superposed to this dynamics, there is a collision process which exchange 
velocities of particles in the same site in a way that momentum is conserved. 

Under diffusive time scaling, assuming local equilibrium, it is not difficult to 
show that the evolution of the conserved quantities is described by the parabolic 
equations 

^ vev 

dtp J +N''Y.VjV VF, {p, p) = Ap, , 

where p stands for the density and p = (pi, . . . ,pd) for the momentum. Fq, . . . , Fd 
are thermodynamical quantities determined by the ergodic properties of the dy- 
namics. 

Consider an initial profile given by {p,p) = (a, I3)+N^^{lpo, (f), where (a, (3) are 
appropriate constants. Expanding the solution of the previous equations around 
(a, /3) and assuming that the first component ipo does not depend on space, we 
obtain that the momentum should evolve according to the incompressible Navier- 
Stokes equation 

div if — , 

dtife = Aadiip'^ + Aiip-Vipi + A2di\ip\'^ + Aipi , 

for 1 < ^ < d, where Aq, Ai, A2 are model-dependent constants. This is the content 
of the main theorem of the article. We prove that under an appropriate time scale 
the normalized empirical measures associated to the momentum converge to the 
solution of the above incompressible Navier-Stokes equation. 

The proof relies on the relative entropy method introduced by Yau [7]. We 
show that the entropy of the state of the process with respect to a slowly varying 
parameter Gibbs state is small in a finite time interval provided the solution of the 
incompressible Navier-Stokes equation is smooth in this interval. 

To obtain such a bound on the entropy, we compute its time derivative which 
can be expressed in terms of currents. A one block estimate, which requires a 
polynomial bound on the spectral gap of the generator of the process, permits to 
express the currents in terms of the empirical density and momenta. The linear 
part of the functions of the density and momenta cancel; while the second order 
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terms can be estimated by the entropy. We obtain in this way a GronwaU inequahty 
for the relative entropy, which in turn give the required bound. 

The article is organized as follows. In Section [2] we establish the notation and 
state the main results of the article. In Sections [3] and IH we examine the incom- 
pressible limit of an asymmetric long range exclusion process. We state in this 
simpler context some ergodic theorems needed in the proof of the incompressible 
limit of the stochastic lattice gas. In Section [5] we prove the main result of the 
article, while in Section [5] we prove a spectral gap, polynomial in the volume, for 
the generator of a stochastic lattice gas restricted to a finite cube and in Section 
[7] we state an equivalence of ensembles for the canonical measures of lattice gas 
models. 



2. Notation and Results 

Denote by = {0, . . . , N — 1}'^ the d-dimensional torus with N'^ points and 
let V C M'' be a finite set of velocities v — {vi, . . . , Vd)- Assume that V is invariant 
under reflexions and permutations of the coordinates: 

{vi,...,Vi_i,-Vi,Vi+i,...,Vd) and (^^.(i) , . . . , '^^(d)) 

belong to V for all 1 < i < d and all permutations a of {!,..., d} provided 
{vi, . . . , Vd) belongs to V. 

On each site of the discrete d-dimensional torus at most one particle for each 
velocity is allowed. A configuration is denoted by r? = {rjx, x e T^} where r]x = 
{r/(a;,w), v <^ V} and r]{x,v) £ {0,1}, a; G T^, t; G V, is the number of particles 

with velocity v at x. The set of particle configurations is Xn — ({0, 1}^)^". 

The dynamics consists of two parts: long range asymmetric random walks with 
exclusion among particles of the same velocity and binary collisions between par- 
ticles of different velocities. The first part of the dynamics corresponds to the 
evolution of a mesoscopic asymmetric simple exclusion process. The jump law and 
the waiting times are chosen so that the rate of jumping from site x to site x + z 
for a particle with velocity v is pn{z, v), where 

PNiz^^) - ]f^{2 + ■^qM{z,v)}i{zeAM} ■ 

In this formula a > is a fixed parameter, M is a function of N to be chosen later, 
Am is the cube {-M, . . . , M^, Aai is given by 

= ^^'1 (2.1) 

for 1 < i, j < d, and qM{z,v) is any bounded non-negative rate such that 

Am 



Md+i ^ 

zSAa 



qMiz,v)zi = Vi 



ior 1 < i < d and M > 1. A possible choice is qM{z, v) = M~^{z ■ v), where u ■ v 
stands for the inner product in W^. Note that particles with velocity v have mean 
displacement M^^N^"-v. 



4 



J. BELTRAN, C. LANDIM 



The generator £^ of the random walk part of the dynamics acts on local func- 
tions / of the configuration space Xm as 

{^^Nfm vix,v)[l^Tj{x + z,v)]pN{z,v)[fir,--^+^'n-fiv)], 

where 

{ri{y, v) ii w = V and z — x, 
ri{x, v) ii w = V and z = y, 
ri{z,w) otherwise. 

The collision part of the dynamics is described as follows. Denote by Q the set 
of all collisions which preserve momentum: 

Q = {(v, w, v' , w') (z V'^ : V + w = v' + w'} . 

Particles of velocities v and w at the same site collide at rate one and produce two 
particles of velocities v' and w' at that site. The generator £^ is therefore 

where the rate p{y, q,i]), q = {v, w, v' , w'), is given by 

P{y,q,v) = v{y,v)fi{y,w)[l-T]{y,v')][l-i]{y,w')] 

and where the configuration rj^''' , q — {vq, fi, W2, ^3), after the collision is defined as 

y.q/ \ _ j viy^ if ^ y and u — Vj for some < j < 3, 

^ 1^ ri(z,u) otherwise, 

where the index of should be understood modulo 4. 

The generator Cn of the stochastic lattice gas we examine in this article is the 
superposition of the exclusion dynamics with the collisions just introduced: 

Note that time has been speeded up diffusively. Let {ri{t) : t > 0} be the Markov 
process with generator and denote by {S^ ■ t > 0} the semigroup associated 
to £Ar. 

For a probability measure 11 on Xn, denote by the measure on the path space 
D{M.+ , Xn) induced by {r]{t) : t > 0} and the initial measure fi. Expectation with 
respect to is denoted by E^. 

2.1. The invariant states. For each configuration ^ e {0,1}^, denote by Io{£,) 
the mass of ^ and by //j(i^), /c = 1, . . . , d, the momentum of ^: 

vev vev 

Set 7(0 := {Iq{^), ■ ■ ■ ,ld{0)- Assume that the set of velocities V is chosen in 
such a way that the unique quantities conserved by the dynamics Cn are mass and 
momentum: J2xeT% 

Two examples of sets of velocities with this property were proposed by Esposito, 
Marra and Yau [2]. In Model I, V = {±ei, . . . , ±6^}, where {ej,j — 1, . . . ,d} 
stands for the canonical basis of R'^. In Model II, d = 3, w is a root of — 6w^ — 1 
and V contains (l,l,w), all reflections of this vector and all permutations of the 
coordinates, performing a total of 24 vectors since w ^ ±1. 
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For each chemical potential A = (Ao, Ai, . . . , Ad) in R''+^, denote by mx the 
measure on {0, 1}^ given by 

mxiO = ^cxp{a.7(0} , 

where Z{X) is a normalizing constant. Notice that nix is a product measure on 
{0, 1}^, i.e., that the variables {^{v) : v €V} are independent under mx- 
Denote by |J,^^ the product measure on ({0, 1}^)^^ with marginals given by 

IJ,x{v ■ V{x, ■) =0 = 

for each ^ in {0,1}^ and x in T^. Notice that {ri{x,v) : x € '^%,v € V} are 
independent variables under ji^. 

For each A in K'^+^ a simple computation shows that /x^ is an invariant state for 
the Markov process with generator ^^v, that the generator Cj^ is symmetric with 
respect to and that £^ has an adjoint -C^'* in which p]\[{z,v) is replaced by 
p*pjiz, v) — pn{—z, v). In particular, if we denote by >C^'*, ^C^'" the symmetric and 
the anti-symmetric part of C^, we have that 

- 7V^d+2 iT^x,x+z,vf){ll) , 

z&Am 



where 



(T.,.+.,„/)(r?) = v{x,v) [1 - rjix + z,v)] [/(r,-'^+^''') - /(»?)] . 



The expectation under the invariant state fj,^ of the mass and momentum are 

given by 

p(A) := E^.MO] = E^-W 

vev 

Pk{\) := EmMO] = Y'^kS^W ■ 

In this formula ^^(A) denotes the expected value of the density of particles with 
velocity v under mx- 

1 + exp |Ao + 2^k=i ^fe) 
Denote by (p,p)(A) := {p(\),pi{X), . . . ,pd{X)) the map which associates the 
chemical potential to the vector of density and momentum. Note that (p, p) is the 
gradient of the strictly convex function log.Z(A). In particular, {p,p) is one to one. 
In fact, it is possible to prove that {p,p) is a diffeomorphism onto 21 C the 
interior of the convex envelope of {1(0, ^ € {0, 1}^}. Denote by A = (Aq, . . . , A^) : 
21 the inverse of {p,p)- This correspondence permits to parameterize the 

invariant states by the density and the momentum: for each (p, p) in 21 we have a 
product measure = Ma(p,p) °^ l}^)^*^- 
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2.2. Spectral gap. For L > 1 and a configuration t], let I^{x) ~ {Iq{x),..., 
I^{x)) be the average of tlie conserved quantities in a cube of length L centered at 
x: 

I\x) I\x,i^) = ^ E ^(^^)- (2.3) 

Let QJl be the set of all possible values of /^(O) when 77 runs over ({0, 1}^)"^^. 
Obviously QJl is a finite subset of the convex envelope of {![£,) : ^ G {0, 1}^}. The 
set of configurations ({0, 1}^) ^ splits in invariant subsets: For each i in 2Jl, let 

ndi) ■■= {r;e ({0,in^^ : J^(0)=i}. 

For each i in 2Jl, define the canonical measure '•^^ the uniform probability 

measure on 7iL(i). 

Denote by the generator C n restricted to the cube Am without acceleration. 
More precisely, on the state space ({0, 1}^) consider the generator £aj„ — + 
which acts on local functions / : ({0, l}''^)^"^ i~> R as 

= E E ^{^.^)[^'v{y,v)]pN{v-x,v)[f{r^-^yn-f{ri)]. 

\x-y\<M 

yeAj/ ges 

Since the only conserved quantities are the total mass and momentum, the pro- 
cess restricted to each component 7iA/(*) is ergodic. It has therefore a finite spectral 
gap: For each i in QJm, there exists a finite constant C(Af, i) such that 

< C(M,i)(/,(-/:A.,/))..,,, 

for all functions / in L'^{v^^,j^i). Here and below {f;f)i, stands for the variance of 
/ with respect to a measure v and (•, for the scalar product in L'^{v). 

We shall assume that the inverse of the spectral gap increases polynomially in 
the length of the cube: There exists Co > and k > such that 

max C(M,i) < CqM". (2.4) 
We prove this hypothesis in Section [6] for Model I. 

2.3. Incompressible limit. For fc = 0, . . . , d, denote by tt*^'^ the empirical mea- 
sure associated to the fc-th conserved quantity: 

x&T% 

where 6u stands for the Dirac measure concentrated on u. 

Denote by {tt'''^,H) the integral of a test function H with respect to an em- 
pirical measure tt'''^ . To compute £n{tt'''^ ,H), note that Cyik{T]x) vanishes for 

= 0, . . . , d because the collision operators preserve local mass and momentum. In 
particular, Cn{tt''^^,H) = N'^C'j^ {n'''^ , H). To compute /:^{7r'='^, iJ), consider 
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separately the symmetric and the anti-symmetrie part of After two summa- 
tions by parts and a Taylor expansion, we obtain that 



M N'i 



i=i xeTl 



for every smooth function H. In this formula, A stands for the Laplacian. Tx stands 
for the translation by x on the state space Xn so that {Txr]){y, v) — rj{x + y, v) for 
all a;, y in Z'', v in V, and Wj^ , j — 1, . . . ,d, is the current given by 



= liJ^E E qM{z,v)z,v{0,v){l-v{z,v)}- (2.5) 



In the same way, for I < k < d, a long but simple computation shows that 
N'^C^-'i-K'"'^ ,H) = {tt'^^^^AH) + 0{M/N), 

J = l x£Tf, 

where is the current defined by 

Wk'!, = ]^E"^- E qMiz,v)z,rjiO,v){l-7j{z,v)}. (2.6) 

vev zeAm 

The explicit formulas for Cn{t^^'^ iH) permit to predict the hydrodynamic be- 
havior of the system under diffusive scaling assuming local equilibrium. By (12.11) . 
the expectation of the currents Wj"' , W^'^^ under the invariant state /x^ are given 

by 

In this formula and below, x(a) = a(l — a). In view of the previous computation, 
if the conservation of local equilibrium holds, the limiting equation in the diffusive 
regime is expected to be 



M 



M 

vev 

where VF stands for the gradient of F. 

We turn now to the incompressible limit. Note that 9y(0) = 1/2 for all v and 
that 

p(0) = M =: ao, Pk{0) = 1E^'= = 0' 

vev 



8 



J. BELTRAN, C. LANDIM 



where the last identity follows from the symmetry assumptions made on V. There- 
fore, A(ao, 0) = and by Taylor expansion, 

x(6'„(A(ao + e(/3o,e¥'))) = 

x(l/2) - e2|_^a,Ao(ao,0)(^^ + - ^ ^ «fc 9,Afc(ao, 0)(^,} + 0(6^) 

i=0 k=l 1=0 

because x'(l/2) = 0, do0viO) ^ (1/4), dkOy{0) ^ (l/4)wfe, I < k < d. Here de 
stands for the partial derivative with respect to the £-th coordinate. It follows from 
the previous explicit formulas for dk9v{0) that 

a,Ao(ao,0) = A6oj\V\-\ a,Afc(ao,0) = ASkj{j2^l}' 
for 1 < fc < d, < £ < d. In particular, x(6'«(A(ao + efo, ecp))) is equal to 

Due to the symmetry properties of V, 

^WfcWj = B6k,j (2.8) 

where i? = X^^ev ^i- "^^^ denominator in the expression inside braces is thus equal 
to B. 

To investigate the incompressible limit around (ao, 0), fix 6 > and assume that 
a solution of (|2.7p has the form p{t,u) = ao + N^^ifio{t,u), pk{t,u) — N^^ifk{t,u). 
Then, to obtain a non-trivial limit we need to set M = N^^°-^^ to obtain that 
((/3o, is the solution of 

l<l<d. 

To recover the Navier-Stokes equation, we need to introduce some notation re- 
lated to the velocity space V. Let Rk,i,m,n{v) = VkVeVmVn- By the symmetry 
properties of V, if A; ^ £, we have that 

Rk,e,m,n{v) — {Sm,kSn,e + 37n,e6n,k} C , (2-10) 

where C = X^^ev ^fc ~ X^ev ^2- On the other hand, 

Rk,k,m,n{v) — 5m,k5n,k D -f 5m,n{^ — 5m,k} C , (2-11) 

where D = ^^gy ^fe = E«gv ■ 

Assume that (^o(0, w) is constant and that div (^(0, •) — 0. Since dt^^ and A(/3o 
vanish and since V is invariant by reflexion around the origin, replacing v by — w, 
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the first equation in ()2.9p can be rewritten as 

d 

^ ^ Wfc w • V<^fc = . 

fc=l v£V 

By (|2.8|) . this equation becomes 

div = . 

The same argument permits to rewrite the second equations in (|2.9p as 

d 2 

The first term on the right hand side of this expression is equal to 

d 

k,rn.n—l vGV 

It follows from (|2.10p . (|2.1ip and elementary algebra that this expression is equal 
to B^^ times 

{D~3C)diipj + 2Cip-yipi + Cdi\ipe\^ 
because div (p = 0. We recover in this way Navier-Stokes equation 

div (f ~ , 



(2.12) 



where ^ {D - 3C)/B^, Ai = 2C/B-^ and A2 = C/B^. For Model I we get 
Ao = I, Ai = A2 = 0, while for Model II, B = 16+8w^, C = 8+16^^, D = le+Su;^ 
and ^0 vanishes because w is chosen as a root of — Gw^ — 1. 

2.4. Statement of the result. Recall that k stands for the polynomial growth 
rate of the spectral gap. Assume that b < a, 

a + b>l~- , a+{ 6>1--, a + (l + - 6 < 1 . 2.13) 

The first two displayed conditions are needed in the proof of the one-block estimate, 
where the size of the cube cannot be too large. The last condition appears in the 
replacement of expectations with respect to canonical measures by expectations 
with respect to grand canonical measures, where the volume [Am | has to be large. 
It is easy to produce constants a, b > meeting the above requirements. It is 
enough to choose first < a < 1, close enough to 1, and then to find b small 
enough. 

Let ip = ((^1, . . . , (pd) : T'' ^ R'' be a smooth divergence free vector field. Denote 
by (p{t) the solution of (|2.12p with initial condition ip, assumed to be smooth in 
a time interval [0,T]. Denote by the product measure on X^r with chemical 
potential chosen so that 

for 1 < A: < d with (/?o being a constant. This is possible for N large enough since 
is bounded and A(ao,0) =0. 



10 



J. BELTRAN, C. LANDIM 



For two probability measures n, v on Xn, denote by iJAr(/i|z^) the entropy of fi 
with respect to v. 

HN{tJ'\i^) = sup I j fd^i - log j e^dv^ , 

where the supremum is carried over all bounded continuous functions on Xn- We 
are now in a position to state the main theorem of this article. 

Theorem 2.1. Assume conditions (12. 4[) and (|2.13p . Let cp — {cpi, . . . , tpd) ■ T'* 

be a smooth divergence free vector field. Denote by ip{t) the solution of (I2.12[) 
with initial condition ip and assume (p{t, u) to be smooth in [0, T] x T'' for some 
T > 0. Let {fi^ : N > 1} be a sequence of measures on X^q such that H]si{^^\vq) — 
o{N'^~^^). Then, HN{ii^S^\yP) = o(N'^-^^) forO<t<T. 

Corollary 2.2. Under the assumptions of Theorem \2.1l for every < t < T and 
every continuous function F : M., 

^ E F{x/N){La{7^At))^a^} = ^0 / F{u)du, 
J™ ^ E F{x/N)Lk{v^{t)) = / F{u)^k{t.u)du 

in Li(P^iv). 

The corollary is an elementary consequence of the theorem and of the entropy 
inequality. 

3. Mesoscopic asymmetric exclusion processes 

We start with a model with no velocities. The proof is simpler in this context 
and the results stated will needed for the stochastic lattice gas. Denote by 77 the 
configurations of the state space Xf^ = {0, 1}^" so that r\(x) is either or 1 if site 
X is vacant or not. We consider a mesoscopic asymmetric exclusion process on Xjq. 
This is the Markov process whose generator is given by 

{L^f){ri) = Y.^{x){\^^{x + z)\vN{z){f{o-■-^^^)-f{^)\, 

where, 



In this formula a > 0, M, Am are chosen as in the previous section and q(y) = 
sign(y • v) for a fixed vector v £ R"^. On the other hand, a^'^rj is the configuration 
obtained from rj by interchanging the occupation variables r]{x), rj{y): 

{r]{z) iizj^x,y, 
T]{y) if z = a; , 
r]{x) a z = y . 

For a probability measure /i on Ajv, stands for the measure on the path space 
_D(R+, A'at) induced by the Markov process with generator Ljv speeded up by N"^ 
and the initial measure /x. Expectation with respect to is denoted by E^. Denote 
by {S^ ■ t > 0} the semigroup associated to the generator N^Ln- 
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For < a < 1, denote by fi^ the Bernoulli product measure on Xn with density 
a. An elementary computation shows that jj,^ is an invariant state for the Markov 
process with generator Ljv- Moreover, the symmetric and the anti-symmetric part 
of the generator Lff, respectively denoted by L^, LJ^, are given by: 



We investigate in this and in the next section the incompressible limit of this 
model. Consider first the hydrodynamic behavior of the process under diffusive 
scaling. Denote by tt^ the empirical measure associated to a configuration: 

Denote by (tt^, H) the integral of a test function H with respect to an empirical 
measure tt^. To compute N'^Ln{'jt^,H), we consider separately the symmetric 
and the anti-symmetric part of the generator. After two summations by parts and 
a Taylor expansion, we obtain that 

N^L%{n'^,H) = (7r^,Aif) + 0{M/N) . 

On the other hand, after a summation by parts, N^Lff{TT^ ,H) becomes 

ATI-" 1 



where stands for the translation by x on the state space so that {Txr]){y) = 
r]{x + y) for all x, y in T^, and = {Wf^ , W^) is the current given by 

= M^+T E ^(0) {1 - ^(^)} • (3-1) 

zGAm 

The expectation of the current under the invariant state /x^ is 

"(1-")^ E (3-2) 

zEAm 

Since M = N^~'^~^, the limiting equation in the diffusive regime is therefore ex- 
pected to be 

atp + ArS-Vp(l-p) = Ap, 

where 7j = uj q{u) du. 

To investigate the incompressible limit around density 1/2, suppose that a so- 
lution of the previous equation has the form p{t,u) = (1/2) -|- N~^ip(t,u). An 
elementary computation shows that 

dtf = 7 • Vi^^ + A(p . 
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Assume the following conditions on a and 6, which could certainly be relaxed: 
d . f 2i. / 2^ 



d + 2 



<a + b, a + max|2,l + -jfe < 1 , 2(l + -)6<l. (3.3) 

The first assumption, which forbids a large mesoscopic range M, is used in the 
proof of the one-block estimate. The second and third assumptions, which require 
a not too small range M, are used throughout the proof to discard error terms. 

By the same reasons of the previous section, there exist positive constants a, b 
satisfying these assumptions. 

Fix a continuous function ipo : ^ M. Denote by = if{t, u) the solution of 
the nonlinear parabolic equation 

For t > 0, let be the product measure on with marginals given by 

This is possible for N large enough because (p is bounded. Recall that we denote 
by HN{iy\fi) the relative entropy of a probability measure v with respect to fi. 

Theorem 3.1. Assume conditions p.3|) . Fix a smooth function ipo : T"^ M. and 

denote by ip — Lp{t,u) the solution of p.4p with initial condition ipo. Assume ip to 
be smooth in the layer [0, T] x T"^. Let {ji^ '■ N > 1} be a sequence of measures on 
Xn such that HNifi^liya) = o{N'^-'^''). Then, Hjv(/i^S'f ) = o{N''~^'') for all 
0<t<T. 

Fix two bounded functions ipi : T"^ ^ M., i — 1,2, and denote by z^^'* the product 
measures associated to the density profile (1/2) + N^^pi. A second order Taylor 
expansion shows that 

i7Ar(l.^'2|^iV,l) ^ OiN"-^') . 

The assumption on the entropy formulated in the theorem permits therefore to 
distinguish between 7V~''-perturbations of a constant density profile. 

A law of large numbers for the corrected empirical measure follows from this 
result. For a configuration ?/, denote by li'^ {rf) the corrected empirical measure 
defined by 

= n^(77) = ^{,7(x)-l/2}4/Ar 

considered as an element of A^(T''), the space of Radon measures on T'' endowed 
with the weak topology. For i > 0, let Ilf^ = n^(?7t). 



Corollary 3.2. Under the assumptions of Theorem \3.1\. for every < t < T and 
every continuous function F : T'^' ^ R, 

lim (nf,F) = / p{t,u)F{u)du 



The corollary is an elementary consequence of Theorem 13.11 and the entropy 
inequality. 
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4. Incompressible limit of mesoscopic exclusion processes. 

We prove in this section Theorem 13. II Fix a smooth function i^o : T'' ^ M and 
denote by f — ip{t,u) the solution of p.4p with initial condition tpo, supposed to 
be smooth in the time interval [0,T]. Let {//^ : > 1} be a sequence of measures 
on Xpf satisfying the assumptions of Theorem 13.11 

4.1. Entropy, Dirichlet form and Ergodic constants. An elementary com- 
putation shows that the entropy i/jv(A*^ 1/^1/2) '^^ with respect to fJ.1/2 is of 
order Indeed, by the explicit formula for the entropy and by the entropy 

inequality, 




'2 



for all ^ > 0. A Taylor expansion shows that the second term on the right hand 
side is of order jSfd.-'^b, In particular 

N''-'H^{^,''\^i^^,) < Co (4.1) 

for some finite constant Co depending only on 930- 

Let /j^ be the Radon-Nikodym derivative d^^ /dfi^^^ that 

where L*j^ stands for the adjoint of Lj^ in L'^^jL^^^). It follows from (|4.ip and a well 
known estimate on the entropy production (cf. [3], Section V.2) that 

^i/^(/i^5f K/2) + Dn{t^^/,j!')ds < Co (4.2) 

for all > and t > 0. In this formula, Dat stands for the Dirichlet form defined 
as 

where (•, is the scalar product in L'^{ij,^^^). An elementary computation shows 
that 

\x-y\<M 

where 

(v-'^5)(^) = gi^^'^v) - giv) , 

and that D^r is a convex, lower semicontinuous functional. 

Let Lam be the symmetric part of the generator Ljv restricted to the cube Am: 

^LKufm ^ E '?(^)[l-'7(j/)][/(^^^^^) -/('])], (4.3) 

a;,aeAA/ 

\x-y\<M 

and denote by ham.k, Q < K < |Am|, the canonical measure on {0, l}'^" concen- 
trated on the hyperplane with K particles. In the case of the exclusion process. 
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MAaj is just the uniform measure over aU configurations of {0, l}'^*^ with K par- 
ticles. Denote by Dam the Dirichlet form associated to Lam ■ 

\x-y\<M 

where /i stands either for the marginal on Am of the grand canonical measure ^'1^2 
or for a canonical measure ^iAM,K- 

By comparing the Dirichlet form with the Bernoulli-Laplace Dirichlet form, 
in which all jumps are allowed with rate |Aa/|^^ and which is known to have a 
spectral gap of order 1 (cf. we can prove that the spectral gap of Dam is of 
order Af~^. 



4.2. The relative entropy method. The proof of Theorem 13.11 is based on the 
relative entropy method introduced by Yau [7]. Let -0^ — dv^ / dji^^^,- It follows 
from the explicit formulas for the product measure that 



Let HN{t) ^ N'^''-'^Hn{pi-'^ S^\v^) and recaU that we denote the Radon-Niko- 

N 
1/2 



dym derivative dfi^ /d^^,^ by fl^. With the notation just introduced, we have 



that 



HN{t) = N^'-'Jfriog^d^^ 



Theorem 13.11 follows from Gronwall lemma and the following estimate. 

Proposition 4.1. Fix a sequence of measures {/i^ ■ N > 1} satisfying the assump- 
tions of Theorem \3.1\ There exists 7 > such that 

HN{t) < 7 / HN[s)ds + on{1) 



Jo 

for all t < T. 

The proof of Proposition 14.11 is divided in several steps. We begin with a well 
known upper bound for the entropy production (see e.g. [3], Lemma 6.1.4). 
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A long and tedious computation gives that [N"^ L*j^ — dt)i^^ is equal to 

^ 5^(A<^)(t,x/iV)Mx)-l/2} 

Ifi 

+ ]V5^ E E {du.v){t,x/N)[d^^v){t,xlN)T,V,^^{r^) (4.5) 
d 

- ]^ E (ft</')(t,:r/iV){ry(a;) - (1/2) - 7V-V(t, :r/iV)} + o{N''-^') . 

In this formula, W*'^ and stand for 

= 17^ E 9(-^)^,^(0){l-r?(z)}, 



- E -.-.-7(0){l-^(^)}. 

We used the inequalities a > b, a + 2b < 1, which follow from assumptions 
to estimate several terms in the above computation by o{N'^^'^^). The expression 
(1 + eN){du.ip) in the third hne stands for (9„^-(p){l — 4:N^^''ip'^}^^ . Keep in mind 
that En is of order N~'^^. 

If we replace ri{x) - 1/2 in the first term of (|4.5p by ?]{x) - {\/2)~ N^'^Lp{t,x/N) 
and V^^ by V^^ — {\/4)5ij, as N t oo, the expressions added multiphed by 
(1/4)A^^''^'^ converge to 



Therefore, in view of (14. 4|) . (14. 5p . the time derivative of the renormalized entropy 
HN{t) is bounded above by 



4(l + e.)iV- ^ E {du,^)it,x/N)r^Wp'\,,)] (4.6) 



^ E E idu,v){t,x/N)id^^^){t,x/N)T,V,'jir^,) 



E 



^ 5](A^-at^)(i,x/iV){ryt(x)-(l/2)-Ar-V(t,^/iV)}] + OAr(l) , 



where Vj^J{r]) = ^,^^(77) - (1/4)5,,^-. 

We now use the ergodicity to replace the functions W*'^^ and Vj^j by their 
projections on the conserved quantity over mesoscopic cubes. 
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4.3. One block estimate. Recall from the end of Subsection 14.11 that ^iA^f.,K 
stands for the uniform measure over all configurations of {0, 1}^*^ with K particles. 



For 1 < j < d, denote by Fj{K/\AM\) the expected value of the current W*' with 
respect to ^Km-K- An elementary computation shows that 

provided P = K/\Am\ and 7f = Af-f'^+D E.^Am ^M^) = 7^" + O(M-i). 

For a positive integer £ > 1, let ri^{x) be the average number of particles in a 
cube of size £ around x: 

^'(^) = E '^(y) • 

I yex+Ae 



.AI 



For M >1, 1 < j < d, let 
Lemma 4.2. For every t > 0, 1 < j < d and continuous function G : T*^ 



F,(V"(0)) 



Hm E„ 



= 0. 



a- ST? 



Proof. By the entropy inequality and Jensen inequality, the expectation appearing 
in the statement of the lemma is bounded above by 



N 



2b 



AN<i 



N' 



2b 



logE^iN 



exp A 



ds J2 Gix/N)TxVjMivs) 



for every j4 > 0. In view of (|4.ip . to prove the lemma it is enough to show that the 
second term vanishes, as t oo, for any ^4 > 0. Since e'^' < + e~^, it is enough 
to estimate the previous expectation without the absolute value. 

By Feynman-Kac formula and by the variational formula for the largest eigen- 
value of an operator, the second term without the absolute value is bounded above 

by 

sup I AG{xlN) / r,y,-Af/dA*f/2 - ^JV (A*f/2 , /) } , (4.7) 



AN'i 



where the supremum is carried over all density functions / with respect to /i(y2- 

Since the measure /i(y2 translation invariant and since Vj_M depends on the 
configuration only through {77(2) : z G Am}, 



(txVj.m) f dfi^^ 



Vj,M (T-xf) dfl^/ 



Vj-AI fx,M dfli/2 ' 



where fx,Ai = E^N^^lT^xflviz) , z e ^m]- 

On the other hand, by convexity of the Dirichlet form and by translation invari- 
ance of A'i/2i ^r any x, y in Am such that |a: — y| < AI, 
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Therefore, summing over x, y in Aa/, |a: — y| < M and in z in T^, we obtain that 

zGT^ x,yeAM \x-y\<M 
\x-y\<M 

for some universal constant Cq. 

Recall the definition of the Dirichlet forms I? at and Dm, Am introduced above. It 
follows from the previous estimates that the expression ()4.7I) is bounded above by 

^ E Bup{^G(x/7V) / F,,M/dA^f/2 - -j^DM,AMi^^l/2J)} (4.8) 

where the supremum is carried over all densities / with respect to the marginal of 
fly 2 on the cube Am- 

In particular, by projecting the density over each hyperplane with a fixed total 
ntimber of particles and recalling the perturbation theorem on the largest eigenvalue 
of a symmetric operator (Theorem 1.1 of Appendix 3 in [5]), in view of (|4.8p . we 
obtain that (|4.7p is less than or equal to 



Ch)At\\G\\l^^{i-LM,AMr%.M,V,,M)^^^ 



for some finite constant C(7) depending only on 7. Here we need the assumption 
that M'^^'^ <C iV^ to be allowed to apply the Rayleigh expansion. 

Since the generator Lm,Am has a spectral gap of order M^^, ((— ^Af.AAj )~^^",M, 
Vj,A/)^N^ is bounded by CqA-P {Vj,MTVj,M) p,^^^, which is less than or equal to Cq 

M'^-'^. Thus, (gJl) is bounded by CoAt||G||^7V-2'^ because M = N^~^-b^ This 
concludes the proof of the lemma. □ 



For 1 < i,j < d, let 

with the same convention that /3 = K/\Am\- Let w^'^v) = - ^*j(^^(0))- The 
arguments of the proof of Lemma 14.21 shows that for every t>0, l<i,j<d and 
continuous function G : T'' — > M, 



lim E,, 



= 0. 



(4.9) 



xeTi 



The arguments are even simpler due to the absence of the factor iV^^ multiplying 
the sum. 
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By Lemma [4.21 and (|4.9|) . integrating in time ()4.6p . we obtain that the entropy 
HN{t) is less than or equal to 



4(1+£^)7V 



2b 



a 



j=i xeTi, 



^ fdsE^.[j2 E idu.^)i^,x/N)id^^ip){s,x/N)F,^,irjfix)) 

° ij = lx6T'?- 



(4.10) 



47^ 



plus an error term of order ojv(l) for every t <T. 

Recall that xi^-) — o^i^ ~ o). Since iV^^ <C M'', we may replace in the previous 
formula f; (r,f(x)) by -jf'xiv^'ix)) and F,,, (j^f (x)) by {x(?7f (x)) - x(l/2)}(5,,,. 
Moreover, since -^"'^ X^^^gT"^ {9ujif){t,x/N) is of order iV~^ and since 6 < 1/2, we 

may further replace xiVt^ (x)) by xiVt^ (x)) ^x(l/2) in the first term. Finally, since 
for a smooth function G, M^'^J2y:\y-x\<AAG{y/N) - G{x/N)] is of order (M/N)^ 

and since A'PN^~'^ vanishes as t oo, we may replace ris{x) by rif^{x) in the third 
term. After all these replacements and since x(&) — x(l/2) = —[b — (1/2)]^, (|4.10p 
becomes 

4(l + e^)7V^'^ 



2b rt d 

- / d.E^«[E E lf'{du^v){s,x/N){rjf{x)-l/2y 



a 

dsE^. [E E (a„.¥')'(5,^/^){^f (2:) - 1/2}' 



J=l xGTl 



(4.11) 



dsl 



[ E (A^ - 9.^)(s,a;/iV){,7f (x) - (1/2) - N-'cp{s,x/N)} 



The second line of the previous formula is easy to estimate. One can argue that 
it is negative or one can add N~^(p{s, x/N) inside the braces and apply Lemma 
below. The first term in (|4.1ip without the factor (1 + e^) can be written as 



AN' 



2b 



Nd 



+ 



Nd 



d.E^«[E E lf{du,^Ks,x/N){Tjf{x)-l/2-N-''^{s,x/N)} 
>[E E l!'idu,v'Ks,x/N){^r{^)-y2-N-''^{s,x/N)} 



ds] 



4 

Nd 



fdsJ2 E l!HduMi^,x/N)^^{s,x/N) 
Jo ._i 



j=l x£T 



As 1 00, for each fixed j, s, the last term of this expression converges to 
4 Jr^a dujj {duj^){s, u) ip'^{s, u) ~ 0. By Lemma H31 below. the first term is bounded 
by 7o Jq dsHN{s) + ojv(l) for some finite constant 70. In the second term, since 
enN^ vanishes as t 00 and since, by (|3.3p . N^ M, we may replace (1 +£Af)7^^ 
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by 7j. The resulting expression cancels with the third term of ()4.1ip because ip is 
the solution of ()3.4|) . This proves Proposition 14. II and therefore Theorem 13. II 

We conclude this section with an estimate on the variance of the density in terms 
of the relative entropy. 



Lemma 4.3. There exists 70 > such that 

^ J2 {vt'(x)-{l/2)-N->'^{t,x/N)} \ < joHNit) + o^(l) 



forO<t<T. 

Proof. By the entropy inequality the expectation appearing in the statement of the 
lemma is bounded above by 

1 „ _ iV^b r r _ / . N 2, 



for every 7 > 0. By Holder inequality, the second term is less than or equal to 



^ ^ logE,«[exp{7|AM|(ryf(a;)-(l/2)-7V-V(t,:r/iV)) } 



The above expectation is bounded uniformly in N provided 7 is small enough. 
The expression is thus bounded by j~^N'^''M~'^, which concludes the proof of the 
lemma. □ 

5. Proof of the Incompressible limit. 

Fix the reference measure = ly^^^ g-j . Consider a sequence of probability mea- 
sures {fj,^ ■ N > 1} satisfying the assumptions of Theorem 12.11 A straightforward 
argument, similar to the one which led to (|4.ip . shows that 

N'''-''HNi^,^\,y^) < Co 

for some finite constant depending only on ao,(po, tp. 

Denote by fl^ the Radon-Nikodym derivative d^^ / dv^ and recall that 
solves the equation 



where £^ stands for the adjoint of Cm in L'^{v^). By the previous estimate on the 



relative entropy of n with respect to , we get that 

+ f dsD^if^) < CoN''-^\ (5.1) 



where stands for the Dirichlet form: D^if) = (/^^^, {—Cn)!^^^)^^. 

Let = dv^ /dv^ . It follows from the explicit formulas for the product mea- 
sures that 

log^f = y: Ht,x)-iM y: i°g ^^zmf ' 

where X{t, x) := A(ao + A^^Vo, N-^ip{t, x/N)) and 

Z{\) - exp{A-/(0} . 

?e{o,i}^ 
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Let HN{t) = N'^''-'^HNitJ,^Sl'\iy^). With the notation just introduced, we have 
that 

HN{t) = N''-'lfriog^d.r 

Theorem 12.11 foUows from GronwaU lemma and the foUowing estimate. 



Proposition 5.1. Fix a sequence of measures {fi^ '■ N > 1} satisfying the assump- 
tions of Theorem \2.1\ There exists 7 > such that 



ft 

HN{t) < 7 / HNis)ds + ow(l) 



for allQ<t<T. 



The proof of Proposition 15.11 is divided in several steps. We begin with a well 
known upper bound for the entropy production. 

j^H^it) < N^'-'j /f . (5.2) 

Next result is needed in to discard irrelevant terms on the right hand side of the 
previous expression. 

Lemma 5.2. Let G : T'' ^ R a continuous function and {ji^ : N >1} a sequence 
of measures satisfying the assumptions of Theorem \2.1\ Then, 



TV-'* ^ G{x/N)h{r]x{t))\ < HN{t) + 0{N-'') 



for 0<t<T,l<k<d. The lemma remains in force for k — if we replace 
hiVxit)) by h{rjx{t)) ~ oq. 

Proof. Fix 1 < fc < d. We may replace hiVxit)) by hiVxit)) - N-^ipk{t, x/N) 
paying a price of order N^^. It remains to apply the entropy inequality with respect 
to measure vt, which is product, and perform a second order Taylor expansion. □ 

A long and tedious computation gives that {ip^)^-^{C'^ — dt)i!^ is equal to 



M 



■^Y. E {^Vk){t,x/N)h{vx) + ;^ E E idu,^k)(t,x/N)TxWl^ 

k=lx£T% j,k=lxGT% 
d d 

E E E (5.3) 



ij = lA:,£=la;eTl 



^(9,^.)(t,x/7V){4(^.)-^^^^^} + R^{t) + o{N'^-'') 



BN'' 

k=l 

In this formula, Wt'f'' and V/'f''^^ stand for 

^kf = ^^tE^*^ E lM{~z,v)zjr]{0,v){l-r]{z,v)} 

^tt^' = -^^Y^'^''^ E z,Zjr){0,v){l-r]{z,v)} . 
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Since the density is a function of the conserved quantities J, the colhsion part of 
the generator is irrelevant in the previous computation. We used repeatedly Lemma 
15 .21 and the fact that b < a, which follows from (|2.13p . to discard superfluous terms. 
The remainder o{N'^~'^^) should be understood as an expression whose expectation 
with respect to fJ,'^ integrated in time is of order o(A^''~^''), while RN{t) is an 
expression which multiplied by N'^^~'^ is bounded by HN{t) + 0{N^^) in virtue of 
Lemma 15.21 

If we replace IkiVx) in the first term of (|5.3p by IkiVx) ~ fk{t,x/N)N^^ and 
i ^-ti ~ {B/4:)dij6k,e, as t cxo, the expressions added when multiphed 



by {B/4)N' 



2b-d 



converge to 



d 

E 

fc=i 



d 

E 



Therefore, in view of ()5.2|) . (15. 3|) and Lemma 15.21 the time derivative of the 



renormalized entropy Hi\}{t) is bounded above by 



HN{t) 



BNd ^ E 



E 



(5.4) 



Ifi 

-WmY. E E(^«.^^)(^'^/^)(^«.^^)(*'^/^)^-^M^''*'W 



r 47V* 



d 

BN"^ E E 



{^^k - dtVk){t,x/N)[^h{Ti^{t)) 



}] + OAr(l) , 



Where V;^'^^ = V^f^'' - (i?/4)<5,,4,.. 

We now use the ergodicity to replace the functions W^''^ and V^j^'^ by their 

projections on the conserved quantities. For s > and x in Z'^, denote by l'^''is,x) 
the average at time s of the conserved quantities over a cube Am centered at x: 

1 



I^'{s,x) 



|A 



M 



E 



To keep notation simple, let 1^ :— I (s,0). 

Recall the definition of the canonical measures vkmA presented in Subsection[ 
Since we assumed in (|2.4p the global dynamics restricted to a cube of length M to 
have a spectral gap of order M'^ and since a + & > 1 — [2/{d + k)], a + (k — 2/k)& > 
1 — (2/k), repeating the arguments presented in the proof of Lemma 14.21 and taking 
advantage of the estimate (15. ip we derive the so-called one block estimate. In 
this lemma, the collision part of the dynamics, also speeded up by iV^, plays an 
important role. 

Lemma 5.3. For every t > every 1 < j, k < d and every continuous function 
G : [0, T] X T'' ^ R, 



lim sup 

Af— »oo 



I '^'l^ E Gis,x/N)r.{w:f{s)^E.^^^^^^,lW:f]} 

xeTi, 



= 
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Since i^A„,i is the counting measure, 



,M 



= -^VkVjE^^^^ . 77(0, -77(61,1;)] 



vev 



because AmM (^^+1) X^zgAm lM{z,v)zj = Vj. In the previous formula, site ei can 
be replaced by any site of Am different from the origin. Since N'^'' ^ M'^, by 
the equivalence of ensembles, stated in Proposition 17. II below, we can replace the 
expectation with respect to the canonical measure by the expectation with respect 
to the grand canonical measure paying a price of order ojv(l). 
For l<j,k<d, let 

where 9y{-) is defined in (|2.2p . Up to this point, we replaced the first expectation 
in (113) by 



E J2idu,^k){t,x/N)R,,k{I^'{t,x))] + on{1) 



j,fe=i xeTt 



Since A(ao,0) = and d^{0) = 1/2, i?j_fe(ao,0) = -(B/4)(5j_fe. On the other 
hand, since (p is divergence free, 

E E {du,Vk){t,x/N)R,,k{ao,0) = -tt^ E E a:/7V) 



vanishes for each fixed N . We may therefore add this expression to the previous 
expectation to obtain that the first term in (|5.4p is equal to 



4A^^ 



E E (5„,^fe)(i,x/7V){i?,-fe(J*^(i,x))-i?,-fc(ao,0)} 



j,k=i xeTL 



(5.5) 

The same arguments show that we can replace V^^ j ' in the second term of 
by its expectation with respect to the grand canonical measure. The proof 
is even simpler due to the absence of the factor TV^'' in front of the sum. Since 
Rj,kiao,0) = -(1/4)5,- feB, 



E.. 



k,e,Ml 



Si^jlRkAP^P) - Rk.eiao,0) 



^A(p,p) L «J J 

The one-block estimate permits therefore to replace the second expectation in 

by 

d 
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E E idu,^kM.x/N){RkAI^'it,x))-RkAao,0)}\ , (5.6) 

i,k,l=l xeT% 

where {du,(pk,e)i*^^/^) = {duifk){t,x/N){dui'Pe){t,x/N). 

It is now clear that (|5.6p is a term of lower order than (|5.5p . We therefore only 
need to estimate the latter. Fix an arbitrary e > 0. Since Rj^k is a bounded 
function, the integral in (|5.5p when restricted to (t, x) — {ao,0)\ > e is bounded 
above by 
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\C,N-' ^ |/*^(t,x)-(ao,0)|^l , (5.7) 



where Cq is a constant depending on V and (p. In the expression above we may 
replace (ao, 0) by (ag + N~^(pQ, N''^(p{t, x/N)) paying a price of order . Since 
x) belongs to a compact set the expression obtained after replacing is bounded 
above by 



E^]N 



By Lemma 15.41 below, this expression is bounded by ^QHpf(t) + ojv(l) for some 
70 > 0. 

In order to deal with the integral (|5.5p on \I^'^{t,x) — (ao,0)| < e we perform 
a Taylor expansion of Rj^k- The first term in the expansion vanishes because the 
gradient of Rj ^ vanishes at (ap, 0). The contribution of the second order terms is 



4iV 
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Expanding the square, the term in Iq vanishes because (p is divergence free and 
the cross product vanishes because V is symmetric. This sum is therefore equal to 

Yl Yi9n,^k){t,x/N)Y^ Vk Vj Vi W„ l[ . 

Replacing by l[ ~ Lpi{t,x/N)N^^, we may rewrite the previous expression 
as the sum of three kind of terms. The first one, the order term in /, consists 
simply in replacing by LpiN~^ . As N tends to infinity, this term converges to 

4.{D - 3C) ^ /■ 2 4C' 



3 = 1- " 3,k=l- 

An integration by parts shows that this expression vanishes because is divergence 
free. The linear term in / cancels with the last term of ()5.4p because ip is the solution 
of the Navier-Stokes equation (|2.12p . Remains the quadratic term in /, equal to 

j,k=i xeTi, 



because p is divergence free. By Lemma 15.41 below, this expression is bounded by 
7offAr(t) + on{1) for some 70 > 0. 

Finally, we consider the remainder in the Taylor expansion. Since Rj^k is smooth, 
we can choose e small enough for the third derivative of Rj^k to be bounded in an e- 
neighborhood of (ao, 0) by a finite constant Co depending on V and (p. In particular. 
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the remainder is bounded above by 



CqN^ 



E 



l''it,x)-iao,0)\ 



The same arguments used to estimate (|5.7p prove that this expression is bounded 
by 7o-ffAr(t) + ojv(l) for some 70 > 0. This concludes the proof of Proposition l5.ll 

We conclude the section with an estimate repeatedly used in the proof of Propo- 
sition [STlJ We need here again the assumption that N'^^ ^ Af*. 



Lemma 5.4. There exists 70 > such that 



Y] < 7oi^Jv(^) + on{1) 



for l<j<d, 0<t<T. The statement remains in force for j — if 
(pj{t, x/N)N^^ is replaced by uq + lpqN^''. 

Proof. By the entropy inequality the expectation appearing in the statement of the 
lemma is bounded above by 



1 

7 



HNit) 



2b 



logE . 



exp{7 J2 (^f (^'^) 



iPj{t,x/N)\^ 



xeTi 



)} 



for every 7 > 0. By Holder inequality, the second term is less than or equal to 



1 E logE-f 



exp 



{7|AMl{/f«,.)-«2}' 



The above expectation is bounded uniformly in N provided 7 is small enough. 
The expression is thus bounded by 'y~^N'^''AI^'^, which concludes the proof of the 
lemma. □ 



6. Spectral gap for stochastic lattice gases 

We prove in this section a spectral gap of polynomial order for the generator 
of the stochastic lattice gas. We consider a slightly different process, in which the 
exclusion dynamics allows particles to jump to any site of Am at rate M~^'^'^^\ We 
do not require, therefore, the jump to be of size smaller than M. Of course, the 
Dirichlet forms of both dynamics are equivalent and the result stated in Proposition 
16.11 extends to the original dynamics. 

Fix M > 1 and consider the process restricted to the cube Am without the factor 
N^. The generator of the process, denoted by Cm, can be written as Cfj + £%j, 
where 

i^MfM = jJ^Y. E Vi^,v)[l-r^iz,v)][f{r,^'^'n-fiv)], 
and p{y, q, rj) is defined at the beginning of Section [2] 
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For each fixed i in QJm , recall that we denote by i^Am A the invariant measure con- 
centrated on configurations r] of ({0, 1}^)^" such that I {rj) = i. An elementary 
computations shows that 

= E E....[{m'''n-mr], (6.1) 

qeQxGAM 

Denote by E^[f; f] the variance of / with respect to a measure v and by (•, •)!^ 
the inner product in L'^{y). 

Proposition 6.1. There exists a finite constant C2, depending only on V, Q, such 
that 

for all f in L^(i/AM,i), all i in 2Jm and all M > 1. 

The proof of this proposition relies on estimates on the Dirichlet forms associated 
to C^^ and Cl^. Denote by the generator of a dynamics in which collisions 
between particles at different sites are allowed: 

= -i^ E E p(^^ 1' ^) {/(^"'') - /(^)} ' 

xi ,...,XiSAM gG Q 

where, for q = {u, v, u', v'), x — {xi, . . . , X4), 

P{x,q,r]) = j^^^{u)-q^^{v)[l-ria:^{u')\[l-r]^^{v')] 

and where jy*"'^ is the configuration r) in which the occupation variables 'qxliu), 
ilx2{v), Vx^iu'), Vxiiv') arc flipped. 

The first lemma of this section states that the Dirichlet forms associated to 
Cj^ and to are comparable and that the Dirichlet form of the conditional 
expectation of a function with respect to the total number of particles with fixed 
velocity can be estimated by the Dirichlet form of the original function. For each 
u in V, let Ky be the total number of particles with velocity v in Am: 

Ky = Ky{r]) = E Vi^^''") ■ 
Lemma 6.2. There exists a finite constant C2, depending only on V, such that 

for every f in L'^{uAM,i)- Moreover, let 
Then, 

{-CmF,F)v^^ ^ < (-/^M/,/),/A„,i 
for every f in L'^{vAM,i)- 

Proof. An elementary computation shows that 

qEQxi,...,X4QAM 
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for / in L'^{h'\^j ,i). Fix q and x. We construct a path from ^ to ^"''^ with jumps 
and coUisions of particles in the same site in the following way. Assume that the 
set V has been ordered: V = {vi, . . . ,Vn} and, without loss of generality, that 
q = (ui, . . . ,W4). We first exchange the occupation variable S,x2{v2), S,xi{v2)', than 
^.x^ivs), S.xiiv3) and finally ^2:4(^4), ^xiiv^)- At this point we may perform the 
collision at site xi and move back the particles and holes to their final positions in 
the reversed order. 

The total length of the path is at most 7. Denote by Co — ,Ci — ^^'"^ 

the successive configurations. Writing {/(C'"''?) - /(Ol as - fiCj)}, 

applying Schwarz inequality, reversing the order of the summations and estimating 
the total number of configurations whose path jumps from Q to Cj+ii we obtain 
that for each q = (vi, . . . , 114), 



E 



Xi,...,X4£Aa 



pix,q,0{f{C-'')-fim' 



i=l x.yeAM 



Co|Aa/| ^ x^i/A, 
xGAm 



[pix,q,o{fie^'')-fior 



for some finite constant Cq. In particular, summing over q in Q and dividing by 
2|AmP, we obtain that 



< 



C2 



|A 



M 



E E 



qeQxeAM 

for some finite constant C2 depending on V. This expression is bounded by 

in view of (|6.ip . This concludes the proof of the first statement of the lemma. We 
turn now to the second. 

Fix X in Am and q = {u, v, w', v') in Q. An elementary computation shows that 



1 



E 



E 



Xi,...,X4£Am 



p{x,q,m{e''')-f{m I {Kviv) -v^V} 



where 



Zg{r,) = K^{r,)K„{r,)i\AM\~K^'mM\~K,') 



In particular, by Schwarz inequality. 



E,, 



[p{x,q,7j){F{rj^^'')^F{rj)} 
P{x,q,v) 



< 



E 



E, 



Xl ,...,X4eAjv 



E 



p{x, q, C){/(r''') - f{i)Y {KM : « e V} 
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Taking conditional expectation with respect to {if 1,(77) : v G V}, summing over x q 
and dividing by 2, the previous expression becomes 

This concludes the proof of the lemma. □ 

Proof of Proposition 16'. II Fix M > 1, i in 2Jm and a function / in L'^ivf^^j i). 
Denote the conditional expectation of / with respect to K — {K.^^v S V} by 
F{K): 

F{K) - E,,^^,[f\{K,„v^V}\ . 

By orthogonality, 

E.,,,Af:f] - E,..,A{f-F^K)Y] + ^-A,,.[^;^] • (6.2) 

Using only the exclusion part of the dynamics, since particles jump uniformly over 
the cube Am with rate M-^'^+'^\ by [S], 

for some finite universal constant Co- 

To estimate the second piece on the right hand side of ()6.2|) . note that the 
exclusion part is irrelevant, while the Dirichlet form associated to the collision part 
can be written as 

qeQ 

Denote by M,i the state space of velocities on Am, 



M,i 



{fsr=(ifi,...,if,):^K=|AM|i} 



and by vAmA the invariant state VAM.i projected on Uma- An elementary compu- 
tation shows that 

for some renormalizing constant Zma- 

Consider from now on model I. Set IAmI* '■— (lo, ■ ■ ■ , Id)- Suppose without 
loss of generality that K-e < Ke for I < j < d and let Kj = K-e ■ Since 
Ij — Kg — K-e 1 K can be recovered from i, Kj. We may therefore ignore 
{Kei , . . • , Kea) and assume that {Ki, . . . , Kd) is evolving on the hyperplane 

H = HM,i = {(A'l, ...,Kd):K,> 0, 2^ifj = Iq - E ^4 " 

j 3 

On the set H. the measure vamA becomes 
For Q<n< \Km\ let 



|AAf 

'<3+I3 



K{a) = ( |AM|-(« + n) x _ 

^ ^ ^ 1 + a y V l + (a + n) y ^ ' 
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hn is a strictly convex, strictly decreasing function in the interval [0, |Aa/| — n\. 
Moreover, for n < m, hn{a) < hm{a) on the interval [0, |AAf| — m]. 

Denote by dj the configuration of N'' with a unique particle at coordinate j. An 
elementary computation shows that 



< 1 if and only if hj^ {Kj) < hi^ (Kk) , (6.4) 



where summation is understood componentwise. 

Denote by K an ordered solution of (|6.7|) below and fix a function F in L'^{DAM,i)- 
We have that 

E,,^^,AF-F] < E,,^^ ,[{F{K)~F{k)r] . (6.5) 

For each K in the hyperplane Ti., consider the following infinite path. Let K'^ ~ 
K and assume that • • • , have been defined. Let jo, such that 

hi,, « - 1) = mmhi, {Ki - f ) , hj^^ {K^J = max/i,^ (Kf) . (6.6) 

If is a solution of ^J}, (/i/,^ (Xf^ - 1) > (-^jo))' ^^^^ ^'^^^ = -^^5 otherwise, 
let K'^+^ = K'^-^ka +5io- In this latter case, by ([Ell), VK^,,i{K'^) < i?Aj„,i(K'^+^). 
Since Ti. is finite and since i>AM,ii^^) decreases whenever is not a solution of 
(|6.7p , the path reaches eventually a solution. The path can therefore be written as 
{K°, ...,K'^°,K'^°,...), where K'''" solves (gJl) and Da,,4K"') < D^mA^") for 
< m < n < Iq. 

By the end of the proof of Lemma below, there is a path from to K of 
length less than or equal to d/2, passing only by solutions of (|6.7p and such that 
all configurations visited have the same probability. Juxtaposing the two previous 
paths, we obtain the path T{K, K) ^ {K ^ K° , . . . , K*^" , . . . , K^"" = K), where 
^ = £k stands for the total length of the path. By construction, the probability of 
the configurations visited is non decreasing. 

We are now ready to estimate the right hand side of (j6.5p . By Schwarz inequality, 

E,,^^^^^[{F{K)-FiK)r] 

< J2i?A,,4K)iK ^{^^(lf,+i)-F(lf,)}2 

K j=0 

Since we just need a polynomial bound on the spectral gap and since this method 
can not provide a sharp estimate, we bound the length of a path £k by the total 
number of configurations \H\ < |Am|''. On the other hand, since D\j^,,i{K) < 
P/^j^j _i(K j) , we may replace the former by the latter. Finally, inverting the order of 
summations and estimating the total number of configurations which contains in 
its path to K a fixed couple Kj, Kj+i by the total number of configurations, we 
get that the previous expression is less than or equal to 

CoA^^'^E E i^K.ui{K){F{L)~F{K)Y 

K Lr^K 

for some universal constant Cq. In this formula, the second sum is carried over all 
configurations L which can be obtained from K by letting a particle jump from a 
site to another: L = K — X3j -\- for some j ^ k. By the explicit formula for the 
Dirichlet form of F derived above, this expression is less than or equal to 

C,M^'i+^''\F,-Cl,F),._^. 
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It remains to apply Lemma [621 to conclude the proof of the spectral gap. □ 
We conclude this section with a result used in the proof of Proposition 16.11 



Lemma 6.3. Fix i in QJm such that h < ■ ■ ■ < Id, IAmI* '■— {lo, ■ ■ ■ , Id) ■ The 

system of equations 

' hi, {Kk - 1) > hi^ (K,) for all l<k,j <d 
2E^'^=^"-E^^ (6-7) 

j j 

0<K,<\Am\, l<j<d, 

has a solution such that Kd < • • ■ < Ki. Moreover, if K , L are two solutions of 
(|6.7p . then \Kj — Lj\ < 1 for all 1 < j < d and Pj^j^,i i{K) = i?/^j^j i(L). 

Proof. To prove the existence of a solution, recall from (16. 4p that 

''^^t{K^^'^ if and only if hj^{Kj) < hi,{Kk ~ I) . (6.8) 

Consider a configuration K* which maximizes the probability VKM,i- The inequality 
on the left hand side of the previous displayed formula is satisfied for all j, k. In 
particular, K* solves (|6.7p . 

Fix a solution of (|6.7p . We claim that Kj < Ki if li < Ij. Assume by contra- 
diction that Ki < Kj. In this case 

hu{K,) < hi^{K,-l) < hi^(K,-l) < hi^{K,), 

which is a contradiction. Here, the first inequality follows from the first property in 
(|6.7p of K, the second from the fact that hj. < hj. and the last from the relation 
Ki <Kj~l. 

Suppose that li — Ij for some i < j and that Ki < Kj for a solution K of 
Let k be such that Kk = Kk for k ^ i, j; K, = Kj, Kj = Ki. It is 
easy to check that K is also a solution of (|6.7p . This observation together with the 
estimate derived in the previous paragraph show that there exists a solution K of 
with Kd<-- <Ki. 
Finally, let K, L be two solutions of (|6.7p . Suppose by contradiction that 
Lj ^ Kj ^ for some j. Since Kj — Lj, there exists i such that Ki < Li. 
In particular, 

hu{L,-l) < hi^{Ki) < hi^{Kj^l) < hi^{Lj) < hj^{L, - I) . 

The first and third inequalities follow from the fact that hj. , hj . are strictly de- 
creasing functions and the relations Ki < Li, Lj < Kj ~ 1; while the second and 
fourth inequalities follow from the property of K, L. This proves the first property 
of K, L. 

To prove the second property of K, L, consider a path from K to L: K = 
M°, . . . , = X, for each < i < £, M''+^ ^ M' + dj - dk for some j = j{i), 
k = k{i). It is not difficult to show that there exists such a path with I < d/2, Af* 
solving (|6.7p for all i. 

Fix i and let = M, M'+^ ^ M + d^-dk, M\ M'+^ solving (gJl). Since M 
solves (p?7)) . hi. (Mj) < hi, {Mk - 1). Using now that M + dj - dk solves we 
obtain the reverse inequality so that hi^ (Mj) = hi, {Mk — 1). In particular, in view 
of (jO)) . j>A„,i(M*) = DAM,i{M'^^). This concludes the proof of the lemma. □ 
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7. Equivalence of ensembles 

We prove in this section the equivalence of ensembles for the stochastic lattice 
gas introduced in Section [51 Recall the definition of the set QJ^ and of the canonical 
measures vk^a- Notice that for every A in W^'^^ 

For in 2t the expectation of the one site random variable I{r]x) under the 

product measure /J'^pp^ is equal to {p,p). It defines a map from 21 to the set of 

probability measures on ({0,1}^)^^. Since this map is uniformly continuous it 
may be extended continuously to the closure of 21. For each i € QJ^ denote by 
MA(i) ^^"^ corresponding product measure by this map. Hence we have a one to one 
correspondence between the canonical measures {vj^^.i : i € 2Jl} and the so-called 
grand canonical measures {MA(i) ■ * ^ ^i}- 

Let {g; /)^ stand for the covariance of g, / with respect to ji: (g; /)^ = E^^lfg] — 
Ep_[f]Ep_[g] and {f,g)^ for the inner product in L?{ij)- 

Proposition 7.1. Fix a cube C A^. For each i e QJ^ denote by the projection 
of the canonical measure vp^^ i on Ag and by fj,^ the projection of the grand canonical 
measure MA(i) - Then, there exists a finite constant C{£,V), depending only 

on £ and V, such that 

\EAf]~E^4f]\<^jXJ'{f-J)li' 

for every f : ({0,1}^)'^' ^ 

Proof. Since is absolutely continuous with respect to /z^, by Schwarz inequality, 
\EAf]-E,4f]\ = \EA{^){f-E,4f])[ 

dv\ dv^ 1/2 1/2 

- ^ V V '^^^ ■ 

Since MA(i) ^ product measure, for any ^ in ({0, 1}^)'^'^, 

d^ ^ Miw[E.gA.\A.Jfe) = |AL|i-|A,|/^(0] 
dl^'^^ M:^(,)[E.eA,Jfe) = |AL|i] 

where J^(^) = \Ae\~'^ J^xeAe ^(^^)- Under MA(i)' ^iVx) are i.i.d. random variables 
taking a finite number of values. By Theorem VII. 12 in 01 j there exists a finite 
constant Co{£,V), depending only on £ and V, such that, 

dv''' 1 
r{C)-l <Co(^,V). 



uniformly in ^. This concludes the proof of the lemma. □ 
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